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Abstract 

While it is commonly believed that there is a direct transition from 
the hadronic to a quark gluon phase at high temperature, it would be 
prejudicial to rule out a sequence of dynamically generated interme- 
diate scales. Using as guide, an effective lagrangian with unconfined 
gluons and constituent quarks, interacting with a chiral multiplet, we 
examine a scenario in which the system undergoes first-order transi- 
tions at T comp , the compositeness scale of the pions, at T x , the scale 
for spontaneous chiral symmetry breaking, and at T c , the confinement 
temperature. 

We find that at current energies, it is likely that the formation 
temperature of the plasma, To < T comp , and that this is therefore a 
quark gluon pion plasma (QGPP) rather than the usual quark gluon 
plasma (QGP). We propose some dilepton-related signatures of this 
scenario. 
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We know that quarks and gluons are confined as hadrons and chiral sym- 
metry is spontaneously broken. We do not however know if these two phe- 
nomena set in at an identical temperature. It is unlikely that the chiral 
symmetry restoration energy/temperature scale is lower than that of con- 
finement, because if it were, hadrons would show parity doubling below the 
confinement but above the chiral breaking scale. This is not seen either exper- 
imentally or in finite temperature lattice simulations. Determining whether it 
is higher even, in the simplest case of QCD with a two flavour SU 2 (L) x SU2(R) 
chiral symmetry is not straight-forward. Indeed, while there is a bonafide 
order parameter for chiral symmetry breaking — namely, the mass of the 
constituent quark — the Wilson loop in the presence of dynamical quarks, is 
not a valid order parameter for confinement — and there are no other known 
candidates. However, by looking at energy density or specific heat, we can 
get a fair idea of the change in the number of operational degrees of free- 
dom or particle modes with temperature. The relevant lattice calculations 
indicate that the drop from the large number of degrees of freedom in the 
QGP/QGPP phase to a few degrees of freedom in the hadronic one takes 
place in one broad step in temperature, suggesting that the two transitions 
are close but not necessarily identical, i.e. that T x ~ T c . 

In the same spirit, we could ask whether pions necessarily come apart the 
moment chiral symmetry is restored, or equivalently, whether the composite- 
ness scale for the pion, T comp , coincides with T x . The following analogy with 
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superconductivity (SC) leads us to conclude that it may not: 



In the usual BCS theory, the formation of Cooper pairs at T crit coin- 
cides with the appearance of a SC order parameter. This corresponds to 
Tcomp = T crit , as there are no Cooper pairs at scales above T crit . Here, T comp 
is the compositeness scale for Cooper pairs, which are normally described as 
undergoing momentum space pairing, the interelectron distance being much 
less than the size of a Cooper pair. This is indeed what happens for a weak 
pairing interaction. For a much stronger interaction, Cooper pairs may ac- 
tually form by real space pairing. The Cooper pair size is then smaller than 
the interelectron distance and it is possible that, even when SC is lost, i.e. 
T > Tarn, pairs continue to exist simply as bound states. 

Evidence for the related precursor phenomena for the chiral transition had 
been suggested long ago JyJ . Recently, there have been several calculations of 



the QCD counterpart of the non-condensed, paired, pseudogap phase of high 
T c superconductivity fl5H , which suggest that the strong pairing situation is 
likely for the strong interactions. 

We shall accordingly take T comp > T x , and argue that heavy- ion collisions 
are the ideal platform for testing this assumption. In doing so, we shall use 
the following SU2(L) x ST^i?) effective lagrangian, in which the quarks are 
coupled to a chiral multiplet, [n, a] [Q, [| |J, to keep track of particle modes 
at various scales: 
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L = --G%G% - £ i, (D + g y (a + ^ 5 rvf)) V> 
-\{9,af - \{d^f - \^ + vr 2 ) - \(a 2 + Tf 2 ) 2 + const (1) 

This has chiral symmetry-breaking but no confinement. The masses of 
the scalar /pseudoscalars and fermions, obtained by minimizing the potentials 
above, are fi 2 = —A < o > 2 and m 2 = 2A < a > 2 respectively. This theory 
reproduces several broad features of the strong interaction at the mean field 
level ||, indicating that it is, perhpas, a reasonable guide to the physics 
at intermediate scales. More specifically, 



1. It provides a nucleon, which is realized as a soliton with quarks bound 
in a skyrmion configuration || [5]. Such a nucleon: 

(a) Gives a natural explanation for the 'Proton spin puzzle': Quarks 
in the background fields are in a spin, 'isospin' singlet state in 
which the quark spin operator averages to zero. If we collectively 
quantize the soliton to get states of good spin and isospin, the 
quark spin operator picks up a small, non-zero expectation value 

(b) Seems to naturally produce the Gottfried sum rule ||. 

(c) Yields, from first principles (albeit with some drastic QCD evolu- 
tion), a set of structure functions for the nucleon which are close 
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to the experimental ones |Q. 

2. As a finite temperature field theory, it yields screening masses that 
match with those obtained from the lattice simulation of finite temper- 



ature QCD with dynamical quarks [fO 



3. It gives a consistent equation of state for strongly interacting matter 
at all densities [ITTL [|. 



This lagrangian and the above discussion on temperature scales lead us to 
consider the following scenario: 



1. For T > T comp , we have gluons and massless current quarks, in a pure 
QGP phase. 

2. For T comp > T > T x , we have a QGPP phase for which the particle 
mode count is (i) 16 modes from 8 massless gluons, (ii) 24 from 2 fla- 
vors of massless current quarks and anti-quarks, and (iii) 3 degenerate 
pions, which are present, in this temperature range, in the form of non- 
Goldstone-boson bound states, which we assume to be massless. (We 
neglect the a, which would otherwise contribute one more mode). 

The operational degrees of freedom in this phase are then 40 — three 
more than in the QGP phase. This increase may be visible as a slight 
bump on an 'Energy Density' vs 'Temperature' plot, in a high resolution 
lattice simulation. 
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3. For T x > T > T c , chiral symmetry is spontaneously broken and the 
degrees of freedom are (i) gluons, (ii) constituent quarks, which acquire 
a constituent mass of 300 — 400 MeV from the spontaneous breaking of 
X-symmetry and are, as a result, Boltzman suppressed given that T x ~ 
T c is expected [|l| to be O(200MeV), (hi) pions as Goldstone bosons, 
and (iv) the a, which also acquires a mass and may be neglected. 

The effective degrees of freedom are thus 19, at this stage. 

4. For T c > T, quarks are confined, and the only massless modes left are 
the 3 Goldstone bosons/pions. 



We can get an estimate of T comp by examining the behaviour of the above 
Lagrangian as a function of energy. We find that at a certain scale, both the 
wave function renormalization constant, Z, for the 7r and a, and the quartic 
scalar interaction, vanish, leaving us with a Yukawa and a scalar mass term. 
The vanishing of Z, or equivalently, of the kinetic term for the scalar and 
pseudoscalars, means that these are no longer dynamical degrees of freedom, 
i.e. they have ceased to exist as composite entities. Eliminating them, via 
their field equation, leaves behind a four-fermi term which gets weaker with 
increasing energy. While these results, gleaned using perturbative RNG, 
cannot be used all the way out to where Z goes to zero (since the Yukawa 
coupling has a Landau singlarity there), they can perhaps give a reasonable 
estimate for the compositeness scale. This works out to be around 700 — 800 
MeV [12], well above the temperatures attained in the present-day heavy 
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ion colliders. This implies that the initial state produced by these is very 
probably a QGPP. 



To investigate this possibility further, we adapt to our scenario, the con- 
siderations of 0, which gives a very transparent and simple treatment of the 
underlying physics. In particular, following [Q], we treat each of our transi- 
tions as first order, with each producing a mixed phase in which the number 
of degrees of freedom changes continuously, as one phase gives way to the 
other. 

A plasma which thermalizes at a temperature, T , where T comp < T <T X , 
a proper time, to, after nuclear impact, then evolves as follows: 



1. It Bjorken expands and cools [13], reaching T x at time, t x \ = (T /T x ) 3 t . 

2. At T x , it undergoes a first-order transition from the quark-gluon-pion 
to the chirally broken, gluon-pion phase, from t x ± to rit x x, where n = 
40/19 is the ratio of the degrees of freedom in the initial and final 
phases. 

N.B. If we assume that chiral restoration and deconfinement occur si- 
multaneously at T x = T c , the quark phase lasts in the mixture from 
time t\ = t x x to t 2 = rt\ (where r = 37/3), which is substantially 
longer. 

3. Having passed to the chirally-broken phase, the plasma undergoes a 



second Bjorken expansion and thereby cools from T x to T c . Since T x ~ 
T c , this part of the evolution may, however, be neglected. 

4. At T c , we get a mixture of the gluon-pion and the purely hadronic 
confined phases, which lasts from rxt x \ to t x2 = r 2 r\t x \, where r 2 is 
19/3. 

5. Finally the pion phase expands from T c to Tf, the freezeout temperature 

at which pions loose thermal contact with one another. This phase 
begins at time t x2 and ends at tf = (T c /Tf)H X 2- 

The basic changes in this scenario when compared to that of JIJ are that, (i) 
the pionic phase starts much earlier, at the initial temperature, To, and not 
after confinement at T c and (ii) the quark phase is suppressed above T c and 
is thus shortened. 

To see how these changes are reflected in the number of dileptons emitted 
by the plasma, we note that these come from ir + n~ and qq annihilations. The 
cross-section for qq — > l + l~, summed over the spin, flavour and color of all 
quarks is 




where the numerical factor, F q = 20/3, if we include just the u and d quarks 
in the flavor sum. 
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The cross-section for n + n — > l + l likewise is, 
<r v [M\ = F w [M](47r/3)(a/M) 2 v /l - 4(m,/M) 2 (l + 2(m ; /M) 2 ) yjl - A(mjM) 2 (3) 
where the pion form factor is given by, 

F W [M] = mj/((m; - M 2 ) 2 + m 2 rj) + (l/4)mj/((mj - M 2 ) 2 + m 2 ,r 2 ,) (4) 

The masses and decay widths of p and p' resonances are m p = 775 MeV, 
m p/ = 1.6 GeV, T p = 155 MeV, r p , = 260 MeV respectively. 

The rate for producing pairs, per unit space-time, with invariant mass 
squared, M 2 , and transverse energy, E T = + M 2 , where p T is the mo- 
mentum transverse to the beam axis, is (for a — q, 7r), 

dNa = a a [M]M 2 _ ^L\ Et K^{E t IT) = A a [M]E T K (E T /T) 

d 4 xdM 2 dE T 4(2vr) 4 V M 2 ' T 0V T/ ; aL J T 0V T/ ; 

(5) 

While using Bjorken's model, it is expedient to write d A x = d 2 XTdytdt where 
t is the proper time and y, the rapidity, of the fluid element. For central 
collisions of equal mass nucleii d 2 XT = where Ra is the nuclear radius. 
We can further perform the t- integral to get: 



dN 



MT 



i ( 

dijdM^dEr = ^aMM)[SE^TX(G[E t /T ] - G[E T /T C \) 



+E T K [E T /T c }(l/2)(r 1 -l)t 2 xi y (6) 



for the quark channel in the multiple transitions (MT) scenario. The first 
term in braces comes from integrating along the cooling curve t = (T /T) 3 t 



from T to T c and the second, from the coexistence line at T x , which we 
approximate as being equal to T c . The shortening of the quark phase is 
reflected in the presence here of r\ = 40/19, as opposed to r = 37/3 for the 
single transition (ST) at T x = T c case. 

The corresponding rate for 7T + 7r~-annihilations is likewise, 



and Ki[x] are the modified Bessel functions. 

In this case, the first term in braces (absent altogether for the ST case) 
results from the cooling of pions in the QGPP, from T and T c , and the 
second, from their cooling in the hadronic phase, from T c to Tf. The latter 
coincides with the ST result if t x2 — > t 2 . The last term comes from integrating 
over time, the volume fraction, f(t), of the pion phase, on the coexistence 
line T x ~ T c . Since pions are present both above and below each of the 
coexistence lines at T x and T c , f(t) = 1 for both these transitions. In the 
approximation, T x = T c , the integral is then simply, 




where, 



G[x] 



x 3 (8 + x 2 )K 3 [x] 



(8) 




(9) 
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The analogous integral in the ST case is, by contrast, (r — l)rt\/2. 



We note that these differential rates further integrate over to the 
closed-form expressions: 

dN MT ( 

= nRlA q (M)TX{Q(To/M)-\H[M/T } - H[M/T C }) 

+ {MlT c ){T /M)- i K 1 [M/T c }{r 1 - 1) j (10) 

where, 

H[x] = x 2 (8 + x 2 )K [x] + 4x(4 + x^K^x] (11) 

and we have translated ratios of times into ratios of temperatures. Similarly, 
dN^ IT 



7r^Ar(M)T 2 ^ 6(T /M) 4 (H[M/T ] - H[M/T C }) + Q( ri r 2 Y(T /M) 4 



dydM 2 

(H[M/T C ] - H[M/Tf]) + (M/T^To/T^KxiM/T^in^Y - 1] \{ 12) 

Finally, the differential rates, dN^ t T /dydM 2 dE T and dN^/dydM 2 for the 
total number of pairs, are obtained by summing the right hand sides of Eqs. (|6|) 
and (0), and Egs-fllCf) and ( [12 ) respectively. 



We can now compare our scenario, in quantitative detail, with that of 
Ref. IJ. For ease of comparison, we use the same parameter- values: to = 
lfm/c, and nR 2 A = 127/m 2 . 

• In Fig.l, we plot dN™ T / dydM 2 dE Tl dN^ IT /dydM 2 dE T and dN% IT / dydM 2 dE T 
vs M, for the combinations of T and T c used in Ref.[|l| but at values 
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of Et, for which (in the first two cases) the ST pion peaks, coming 
from annihilation via the p and p' resonances, are largely extinguished. 
We note that in the MT scenario, the rate displayed reduces in magni- 
tude as Et increases, but the pion peaks survive and stay well above 
the qq contribution. This results from pions being created when the 
plasma forms, as opposed to when it hadronizes. It is consistent with 
the observation in QXJ] that peak extinction at large Et and the ini- 
tial temperature of the pions are distinctly correlated: the smaller this 
temperature, the smaller the Et at which the peaks get extinguished. 

We note that an enhancement in in the production of large Et dilep- 
tons, owing to pions being available at higher temperatures, will in- 
evitably be accompanied by an enhancement in the production of large 
Pt pions, a feature highlighted by Schucraft (TB| long ago. We can add 



that pions present at the boundary of the initial plasma, will, being 
colourless, have no problem escaping from the interaction region. The 
availability of pions at higher temperatures is, incidentally, supported 
by the observation that the HBT size of the initial pion source is smaller 



than if pions formed only after hadronization [17, 18 



We would intuitively expect the difference between the two scenarios 
to depend (among other factors) on the time spent by the system in 
the QGPP phase. This would in turn be determined by the initial 
temperature, T . In Fig. 2, we accordingly plot the ratio of dN tot / dydM 2 
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for the MT and ST scenarios vs T over the range T x < T < T comp , 
for several values of M. We find that the MT rate is 20 - 30% higher, 
when M is not too far from the p, p' resonances, but drops sharply as 
M increases and the drooping pion form factor begins to cut down the 
pion contribution. 

• We have noted that quarks acquiring constituent masses larger than 
T x and decoupling at the chiral phase transition, reduces the lifetime 
of the quark phase above T c . This, in turn, brings down the qq — >■ 
l + l~ rate for small M and E T . In Fig.l, the difference is visible but 
only when T c = 2A0MeV and T = 250MeV. In the total rate, the 
effect is swamped by the pion contribution for those values of M for 
which the pion form factor is appreciable. However, from Fig. 3, we see 
that for slightly larger values of M, the pion contribution disappears 
and the total differential rate dN tot /dydM 2 reduces to just the quark 
contribution, for both the ST and MT cases. Determining the rate 
in this M-window would then give an indication of whether a chiral 
transition has occurred close to but distinctly above T c . 

• The behaviour of these curves as functions of T can be further un- 
derstood as follows: In the ST scenario, a higher T does not change 
the pion contribution to dilepton production. However, it clearly in- 
creases the time spent by the quarks in the cooling phase and thus 
enhances the quark contribution to dileptons, pushing the crossover 
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to the quark dominated regime to smaller M. For the MT scenario, 
both the quark and pion contrbutions to dileptons are enhanced and so 
the crossover to the quark-dominated regime is not expected to change 
much (though it may still move slightly towards smaller M, due to steep 
fall in the pion form factor). For both scenarios, the cooling term in the 
quark contribution increases in magnitude with T while the Maxwell 
co-existence term (which accounts for the difference between the two 
scenarios) remains unchanged. Thus as T increases, the ST and MT 
quark contributions approach each other. 

• Finally, in Fig.4, we plot the ratio of dN tot / dydM 2 for the MT and 
the No Transition (NoT) scenarios vs T . The latter corresponds to 
there being no phase transition, i.e. to the fireball's cooling in the pure 
hadronic phase. It is seen that this ratio is enormous. 

• We conclude that if the dilepton rates increase dramatically, a phase 
transition has probably occurred. If the pion peaks in dN tot /dydM 2 dE T 
vs M survive to high values of Et, the plasma has a good chance of 
being a QGPP. If the peaks do not survive, but rather reduce to a flat 
line, we are quite likely seeing a quark phase. The actual value of this 
constant rate, for the window in M discussed above, will then decide 
whether the transition to the confined hadronic phase occurred directly 
or through a chirally broken phase. 
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In summary, we have argued that there are reasons to believe that a 
QPG may form in heavy-ion collisions, as a result of several, as opposed to 
a single, phase transition. We have explored this possibility on the basis of 
some simple assumptions. We have found that this scenario pushes to higher 
temperatures the advent of the QGP, but opens up, in return, the exciting 
possibilty of a QGPP, containing pions as non-Goldstone-boson bound states, 
well above T x . We have further investigated experimental signatures based 
on dilepton production, which would help distinguish the single- from the 
multiple-transition scenario. Making this distinction would shed light not 
merely on how the QGP forms but on a range of assumptions which lie at 
the very core of our understanding of strong interaction dynamics. 
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Figure Captions 

Figure 1 Plots of dN t ^ T /dydM 2 dE T vs M for three (T , T c )-combinations. 
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For the values of E T in Fig.l.a and Fig.l.b, the pion peaks in the single 
transition rate are essentially non-existent. Also plotted are the quark con- 
tributions for both the MT (thick line) and ST (thin line) scenarios. These 
can be distinguished only for T = 250MeV and T c = 240MeV. The MT 
quark contribution is lower because of the shortening of the quark phase. 

Figure 2 Plots of the ratio of the dN tot /dydM 2 rate vs T for the MT and 
ST scenarios. For the lower M-values considered in these graphs, the pion 
form factor is appreciable. It seen that the MT rate is distinctly higher than 
the ST rate. 

Figure 3 Plots of the ratio of the dN tot /dydM 2 rates for the MT and ST 
scenarios vs T . The pion contribution is essentially absent for the higher 
M-values chosen. Note that the MT rate is now lower than the ST rate. 
The dashed lines represent the ratio of just the quark contributions. Quarks 
contribute fewer dileptons in the MT scenario, where the quark phase is 
shortened by the occurrence of the chiral transition. 

Figure 4 The ratios of the MT rate to the rate if no transition occurred at 
all. 
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